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Abstract 

We determine the range of parameter space of a Interacting Quintessence Model that best fits the 
recent WMAP measurements of Cosmic Microwave Background temperature anisotropies. We only 
consider cosmological models with zero spatial curvature. We show that if the quintessence scalar 
fields decays into cold dark matter at a rate that brings the ratio of matter to dark energy constant 
at late times, the cosmological parameters required to fit the CMB data are: dark energy density 
iix = 0.43 ± 0.12, baryon fraction = 0.08 ± 0.01, slope of the matter power spectrum at large 
scales Us = 0.98±0.02 and Hubble constant Hq = 56±4 km/s/Mpc. The data prefers a dark energy 
component with a dimensionless decay rate parameter = 0.005 and non-interacting models are 
consistent with the data only at the 99.9% confidence level. Using the Bayesian Information 
Criteria we show that this extra parameter fits the data better than models with no interaction. 
The quintessence equation of state parameter is less constrained; i.e., the data sets an upper limit 
Wx < —0.86 at the same level of significance. When the WMAP anisotropy data are combined with 
supernovae data, the density parameter of dark energy increases to ilx — 0.68 while A augments 
to 6.3 X 10“^. Models with quintessence decaying into dark matter provide a clean explanation 
for the coincidence problem and are a viable cosmological model, compatible with observations of 
the CMB, with testable predictions. Accurate measurements of baryon fraction and/or of matter 
density independent of the CMB data, would support/disprove these models. 
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I. INTRODUCTION 


The recent observations of higi redshift supernovae [Jj], Cosmic Microwave Background 


n 


Q. 


(CMB) temperature anisotropies ^ and the shape of the matter power spectrum consis- 
tently support the idea that our Universe is currently undergoing an epoch of accelerated 
expansion Currently, the debate is centered on when did the acceleration actually start 
and what agent is driving it. A variety of models based on at least two matter compo¬ 
nents (baryonic and dark) and one dark energy component (with negative pressure) have 


been suggested -see p. The ACDM model, where a vacuum energy density or cosmolog¬ 
ical constant provides the negative pressure, was the earliest and simplest to be analyzed. 


While this m^^is consistent with the observational data (high redshift supernova Q, CMB 


anisotropies galaxy cluster evolution p), at the fundamental level it fails to be con¬ 

vincing: the vacuum energy density falls below the value predicted by any sensible quantum 


held theory by many orders of magnitude |^, and it unavoidably leads to the coincidence 
problem, i.e., “Why are the vacuum and matter energy densities of precisely the same order 
today?” More sophisticated models replace A by a dynamical dark energy either in the 
form of a scalar held (quintessence), tachyon held, phantom held or Chaplygin gas. These 
models ht the observational data but it is doubtful that they solve the coincidence problem 


S 


Recently, it has be en p roposec 


evolve separately [10, |11[ ll^ 


5 De en p roposea tnat oarK mat 


(IQ) models of references |1C 


that dark matter and dark energy are coupled and do not 



In particular the Interacting Quintessence 


, llll llO, llO]; aside from htting rather well the high redshift 
supernovae data, quite naturally solve the coincidence problem by requiring the ratio of 
matter and dark energy densities to be constant at late times. The coupling between matter 


and quintessence is either motivated by high energy particle physics considerations 


3 


or is constructed W requiring the hnal matter to dark energy ratio to be stable against 
perturbations (lO.h,?j|. Since the nature of dark matter and dark energy are unknown there 
are no physical arguments to exclude their interaction. On the contrary, arguments in favor 


of such interaction have been suggested [l^. As a result of the interaction, the matter 
density drops with the scale factor a{t) more slowly than a~^. 

A slower matter density evolution hts the supernovae data as well as the ACDM concor¬ 


dance model does 


The interaction also alters the age of the Universe, the evolution of 
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matter and radiation perturbations and gives rise to a different matter and radiation power 
spectra. All these effects will be used to set constraints on the decay rate of the scalar held 
using cosmological observations. In this paper, we shall further constrain the Chimento et ah 
model by using the recently WMAP measurements of the cosmic microwave background 
temperature anisotropies. As it turns out, a small but non-vanishing interaction between 
dark matter and dark energy is compatible with the WMAP data with the advantage of 
solving the coincidence problem. To some extent, this was already su gges ted in a recent 
analysis that uses the position of the peaks and troughs of the CMB jlSj to constrain a 
general class of interacting models designed not to strictly solve the coincidence problem, 
but to alleviate it Briehy, the outline of the paper is: in Section II we summarize the 
cosmological model, in Section III we derive the equations of dark matter and dark energy 
density perturbations and hnd the range of parameter space that best hts the observations; 
finally, in Section IV we discuss our main results and present our conclusions. 

II. THE INTERACTING QUINTESSENCE MODEL 


The IQ model considered here has been constructed to solve the coincidence problem by 
introducing a coupling between matter and dark energy; their respective energy densities 
don’t evolve independently. In this paper, we shall simplify the Chimento et al. model Q 
in the sense that the IQ will be assumed to decay into cold dark matter (CDM) and not into 

nn 

baryons, as required by the constraints imposed by local gravity measurements |2^ The 
baryon-photon fluid evolves independently of the CDM and quintessence components. Un¬ 
like |l^, we do not include dissipative effects. In [l^ this scaling was considered only during 
matter domination, so the scalar held would evolve independently of the CDM component 
until it started to decay at some early time. These assumptions facilitate the numerical 
work while they preserve its essential features. Specihcally, the quintessence held (denoted 
by a subscript x) decays into pressureless CDM (subscript c) according to 


dpx 

dt 


+ 3Hpc — 3Hc^ (pc + Px) 1 
+ 3(1 -t- Wx)Hpx = -3Hc^ (pc + Px), 


( 1 ) 


where tCa, < 0 is the equation of state parameter of the dark energy and is a small 
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dimensionless constant parameter that measures the intensity of the interaction. Approaches 
similar (but not identical) to ours have been discussed in [lO, lUl llSl lla, ll7| . Eqs. ((T)) were 
not derived assuming some particle physics model for the interaction, where quintessence is 
described as an scalar held with a given potential. We followed a phenomenological approach 
and instead we have required the Interacting Quintessence Model to solve the coincidence 
problem. We have imposed the dark matter to dark energy density interaction to give a 
dark matter to dark energy ratio constant at late times and stable against perturbations. 
As a result, the shape of the scalar held potential is also hxed. 

Eqs. (HD can be solved by Laplace transforming the system. The result is 


Px{a) = —-[3(c^f]c,o - (1 - c^)^x,o)(a^+- + ^x,o(5'-a‘^-- S'+a‘^+)], 

STrGwe// ^2) 

Pc(fl) = o ^°-[3((1 + Wx + c^)ffc,o + ~ — S'+a'^+)], 

oTTUWeff 


where Wg// = + and S± = —3(1+ tna;/2) =f (3/2)r(;e//. The density parameters 

fIc,o and 123,^0 denote the current values of matter and dark energy, respectively. Solutions 
of Eqs. ((H) are plotted in Fig. [U Solid, dashed, dotted and dot-dashed lines correspond to 
flc, f2a;, VLr and fib, respectively. In panel (a) = 0.1 and there is a short period of baryon 

dominance; this does not happen in panel (b) where = 5 x 10“^. 

As detailed in Ref. the interaction dark matter-dark energy brings the ratio r = pdpx 
to a constant, stable value at late times. From Eqs. CD, it is seen that the evolution of the 
aforesaid ratio is 


dr 

dt 


3Hd 


^2 _|_ 



-f 2 j r 1 


The equation dr/dt = 0, has two stationary solutions, namely. 


r± = -Wx/{2d) - 1 ± [w^/(4c'‘) +Wx/dY^‘^ , 


( 3 ) 


which verify r+r_ = 1 (with r+ > r_). As shown in Fig. |21 the ratio evolves from an 
unstable maximum r+ at early times -with dark matter and quintessence energy densities 
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Figure 1: Redshift evolution of different energy densities. Solid, dashed, dotted and dot-dashed 
lines correspond to Oc, and Qb, respectively. In panel (a) cP = 0.1, and in panel (b) 

P = 5 X 10“^. The following parameters were assumed: = 0.25, Qx,o = 0.7, = 0.05, 

11^,0 = 10“®, and Wx = —0.99. 


scaling as to a stable minimum r_ at late times, where both energy densities scale as 
a^-. In Fig. |2l the smaller the coupling constant the larger the ratio of cold dark matter 
to dark energy in the past and the smaller in the future without signihcantly affecting the 
length of the transition period. We are not suggesting that the Universe is already in the 
late time epoch of constant, stable ratio r_. The value of the asymptotic ratio is determined 
by the strength of the interaction and at present this ratio could be still slowly evolving in 
time. 

In terms of a scalar field description, the second equation of is equivalent to 


d?(j) 

dP 




( 4 ) 


where 0 denotes the dark energy field and 14//(0) is the effective potential. The latter is 
given by 


"^eff =+ {pc +px)/{d(j)/dt). (5) 

If r = constant, the potential has two asymptotic limits: V oc during both matter 

domination and the period of accelerated expansion, and V oc with a > 0 well within 
the radiation dominated period. A detailed study has shown that only potentials that are 
themselves power laws -with positive or negative powers- or exponentials of the scalar held 
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Figure 2: Evolution of the ratio r = pdPx from an unstable maximum toward a stable minimum 
(at late times) for different values of d. We took tq = 0.42 as the current value. 


yield energy densities evolving as power laws of the scale factor [2]| . Potentials with expo¬ 
nential and power-law behavior have been considered extensively in the high energy physics 
literature. Exponential potentials arise as a consequence of Kaluza-Klein type compacti- 
hcations of string theory and in iV = 2 supergravity while inverse power law models arise 
in SUSY QCD (see 1^|2^). Potentials showing both asymptotic behaviors have also been 
studied j2^, but at present there are not satisfactory particle physics model to justify the 
shapes of potentials of this type (24 1. 

From the evolution of the background energy densities it is possible to constrain the 
amplitude of the IQ and CDM coupling. Since Wg// > 0, is conhned to the interval 
0 < < |tCa;|/4. Negative values of d would correspond to a transfer of energy from 

the matter to the quintessence held and might violate the second law of thermodynamics. 
Further constraints can be derived by imposing stability of the interaction to hrst order loop 
corrections 0 . In Fig. El we used the supernova data of Riess et al. Q| to constrain model 
parameters. In the hgure we plot the 68%, 95% and 99.9% conhdence levels of a cosmological 
model after marginalizing over, Wx and the absolute magnitude of SNIa. We set a prior: 
— 1.0 < Wx < —0.6. Variations of the baryon density produce no signihcant differences and 
the Hubble constant is unconstrained by this Hubble test since the absolute luminosity of 
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Figure 3: Joint confidence intervals at 68%j_^95% and 99.9% confidence level of IQM fitted to the 


“gold” sample of SNIa data of Riess et al. 


Q 


Type la supernovae is not accurately measured. The contours are rather parallel to the 
(? axis, i.e., the low redshift evolution of interacting models is not very different from the 
non-interacting ones. 


III. OBSERVATIONAL CONSTRAINTS ON THE MATTER-QUINTESSENCE 
COUPLING 

Primordial nucleosynthesis and Cosmic Microwave Background temperature anisotropies 
provide the best available tools to constrain the physics of the early Universe. By assumption, 
the scalar held decays into dark matter and not into baryons. Since dark matter and 
quintessence density perturbations are coupled to baryon and photons only through gravity, 
there is no transfer of energy or momentum from the scalar held to baryons or radiation. 
The evolution of density perturbations of dark matter and dark energy can be simply derived 
from the energy conservation equation. In the equations below we shall use the conventions 
of j^. In the synchronous gauge. 
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( 6 ) 


6 , = + 

2 a \r 

Oc = 0 , 


while the evolution of dark energy density perturbations is given by 

2 ^ Cl 

- 9 {cl^-w^){l + w^)e^k~‘^ + ?,-c^{5^ + r5c) 

\a / ’ a 

dx = —(1 — 3c^ )-6'3; + ’ 5x — 3-t—- {^ + r)6x, 

’ a 1 + Wx al + Wx 


( 7 ) 


where 6 and 6 denote the density contrast and the divergence of the peculiar velocity held of 
each component, respectively, derivatives are with respect to conformal time and is the 
quintessence sound speed, taken to be unity as for a scalar held with a canonical Lagrangian. 
The interaction introduces the terms with a factor on the right hand side of Eqs. and 
©• Combining these equations, the evolution of density perturbations in the IQ held are 
described by a driven damped harmonic oscillator, where the driving term is the gravitational 
held After a brief transient period, the evolution is dominated by the inhomogeneous 
solution and is insensitive to the initial amplitude. To hnd the model that best hts the 
WMAP data, we have implemented equations 0 , (0 and ( 17 ) into the CMBFAST code Q. 
We used the likelihood code provided by the WMAP team 29| to determine the quality of 
the ht of every model to the data. Since we are introducing a new parameter, the coupling 
between dark matter and dark energy, the parameter space could become degenerate with 
diherent local maxima representing models that ht the data equally well. For this reason, we 
did not use a Monte Carlo Markov Chain approach Q but we run through a grid of models 
on a six-dimensional parameter space. Grids of models are computationally very expensive. 
To make the computations feasible we reduced the parameter space by introducing prior 
information. We imposed two constraints: (1) all models were within the 90% conhdence 

n 

level of the constraint imposed by Big Bang Nucleosynthesis: 0.017 < < 0.027 and 

(2) in all cosmologies the age of the Universe was chosen to be to > 12 Gyr. With these 
requirements, we explore the region of parameter space close to the concordance model. 
We have considered only hat models with no reionization, no gravitational waves and no 
running of the spectral index. We considered a 6-dimensional parameter space and assumed 
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Figure 4: Joint confidence intervals at the 68%, 95% and 99.9% level for pairs of parameters after 
marginalizing over the rest. For convenience the axis is represented using a logarithmic scale 
and it has been cut to (? < 10“^, though models with = 0 have been included in the analysis. 
In panels (a), (b) and (c) models were fit to CMB data alone. In panel (d) we included supernovae 
data of Riess et al. Q| . 

our parameters to be uniformly distributed in the following intervals: Hubble constant 
Hq = [46,90]/cm/s/Mpc, baryon fraction Qh = [0.01,0.12], dark energy = [0.1,0.9], 
slope of the matter power spectrum on large scales Us = [0.95,1.04], dark energy equation 
of state Wx = [—1.0, —0.65] and = [0,0.05]. We took 23, 15, 33, 10, 9 linear subdivisions 
and 22 logarithmic subdivisions of the above intervals, respectively. The likelihood was 
computed using the routines made publicly available by the WMAP team. 

In FigEJwe give conhdence intervals for pairs of parameters after marginalizing over the 
rest. Contours represent the 68%, 95% and 99.9% conhdence levels. Models with = 0 
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have been computed and were included in the analysis. The results were undistinguishable 
from those of (? = 10“^ the last point included in the graphs. The WMAP data sets strong 
upper limits on the quintessence decay rate. A non-zero decay rate is clearly favored by 
the data. At (? ~ 10“^ the contours indicate a steep gradient in the direction of growing 
(?. This behavior is associated with the decreasing fraction of CDM at recombination 
with increasing (?. When the interaction rate is large, the Universe goes through a period 
dynamically dominated by baryons (Fig^). The oscillations on the baryon-photon plasma 
induce large anisotropies in the radiation and those models are strongly disfavored by the 
data. In Fig0J the models £t the data more comfortably with lower values of fla, and Hq than 
in the ACDM concordance model. Interacting models have larger dark energy density in the 
past than non-interacting models, achieving the same rate of accelerated expansion today 
with a smaller Our best model also requires larger baryon fraction since the matter 

density is smaller prior to recombination than in the concordance model, therefore dark 
matter potential wells are shallower and a higher baryon fraction is required to reproduce 
the amplitude of the hrst acoustic peak |^. The mean value of the cosmological parameters 
and their corresponding la conhdence intervals are: fla; = 0.43 ± 0.12, fib = 0.08 ± 0.01, 
Hg = 0.98 ± 0.2 and Hq = 56 ± Akm/s/Mpc. The latter number is not very meaningful 
since the probability distribution of Ho is rather skewed. As it can be seen in FigH^, low 
values of Hq are suppressed very fast. As the height of the hrst acoustic peak scales with 
Vtbh?, high values of baryon fraction are speedily suppressed by the WMAP data, which 
translates into an even faster suppression of low values of Hq. With respect the quintessence 
equation of state, as we did not explore models with Wx < —1, we can only set an upper 
limit Wx < —0.86 at the la confidence level. Finally, as we chose a uniform prior on logc, 
the conhdence interval is not symmetric, resulting: = 0.005 ^qqq3. 

Our main result is that models with interaction are preferred over non-interacting models, 
with the remarkable feature that they require very diherent cosmological parameters than 
the concordance model. The range for fla,, Hq, fib and Wx are not directly comparable to those 
found in since the interaction is diherent and we used diherent priors. Their coupled 
model requires higher values of the Hubble constant when the strength of the interaction 
increases, opposite to the behavior found in Fig. H] Non-interacting models are compatible 
with the data only at the 99.9% conhdence level. Our best ht model (c^ = 5 x 10“^, 
fla; = 0.43, fib = 0.08, Ho = bAkm/s/Mpc, Ug = 0.98, w = —0.99) has a = —2log{C) = 
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Figure 5: Radiation Power Spectrum. The solid line is our best fit model (c^ = 5 x 10“^, Q.x = 
0.43, Ofe = 0.08, Hq = 5Akm/s/Mpc, Ug = 0.98,re = —0.99). Dashed line corresponds to the 
ACDM concordance model and dot-dashed line is QCDM with parameters = 0.5, Qb = 0.07, 
Ho = GQkmIsjMpc, w = —0.75 and Ug = 1.02. 


974 while the best fit for a non-interacting model occurs at hla, = 0.5, VLb = 0.07, Hq = 
QOkm/s/Mpc, w = —0.75, Ug = 1.02 and has = 983. The Bayesian Information Criteria 
defined as BIC = + k\ogN 3^, that penalizes the inclusion of additional parameters to 


describe data of small size (in this case the number of independent data points is iV = 
and k, the number of model parameters, is 5 in the non-IQ model, and 6 in the model with 
interaction), gives ABIC = —2 which can be considered as positive evidence in favor of 
including this additional parameter to describe the data. The ACDM concordance model 


was deduced by fitting a different set of parameters to WMAP data |29| and the results 
are not directly comparable to ours. For completeness, let us mention that the concordance 
model with Da = 0.72, = 0.049, Hq = Q8 Km/s/Mpc and Ug = 0.97 has a fit of = 972 

is positively better than ours if the amplitude of the matter power spectrum and the redshift 
of reionization are included as parameters. If only the overall normalization of the power 
spectrum is included, the £t is = 990. In this case, the BIC would give ABIC = —9 that 
must be taken as strong indication that the interaction improves the £t. 

The likelihood curves of Fig. HJ seems to suggest that our model is ruled out by ob- 
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servations since, for example, luminosity distance estimates from high redshift supernovae 


indicate that > 0.6 at the 95% conhdence level [l|. However, analysis of the temperature- 
luminosity relation of distant clusters observed by XMM-Newton and Chandra satellites 


appear to be consistent with hlc up to 0.8 [33|. Moreover, with no priors on hlc and tCa,, (so 
phantom models are included in the analysis) values of the CDM fraction as high as ours 
are consistent with SNIa data 3^. Contours in Fig. ^ are rather parallel to the r 23 ;-axis, 
while at the redshifts probed by supernovae interacting models behave as non-interacting, 
and in the same plane contours are parallel to the c^-axis (see Fig. El). In Fig. mi we plot the 
conhdence intervals combining both WMAP and high redshift supernovae data. Contours 
are shifted to large values of dark energy, the reason being that WMAP data constraints the 
coupling while the supernovae data is insensitive to it. For illustrative purposes, in Fig. 
Elwe compare our best ht model (solid line) with the concordance model (dashed line) and 
the best quintessence model with = 0 (dot-dashed). All models are rather smooth com¬ 
pared with the rigging present in the data with the excess coming from similar regions in 
/-space: / ~ 120, 200, 350. The discrepancy among the models is clearer at the hrst acoustic 
peak. This rigging introduces a high degeneracy among cosmological parameters since very 
different models £t the data with similar per degree of freedom. When adding other data 
sets, like the supernovae in Fig. Et:, the fraction of dark energy increases dramatically. With 
respect to the Hubble constant, the value is less certain since several groups advocate values 


close to QOkm/s/Mpc j^, and even as low as 50km/s/Mpc 

There is a signihcant difference between the concordance model and our best fit model. 
Since the latter requires lower Hubble constant and dark energy density, it generates a smaller 
Integrated Sachs Wolfe effect, responsible for the rise of the radiation power spectrum at 
I < 10. This is a generic feature of this type of IQ models and, therefore, the low amplitude 
of the measured quadrupole and octupole is less of a problem in models with decaying dark 
energy than in the concordance model. 


IV. DISCUSSION 

We have shown that a model where the coincidence problem is solved by the decay of the 
quintessence scalar held into cold dark matter is fully compatible with the WMAP data. The 
best model, ~ 5 x 10“^, hts the data significantly better than models with no interaction. 
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Our best fit model requires cosmological parameters, in particular and Ho, that are very 
different from the concordance model. Our models are highly degenerate in the — plane 
with contours almost parallel to the Oj, axis (see Fig. HD but including prior information 
from supernovae data shifts this parameter close to the concordance model Oj, ~ 0.68. We 
wish to emphasize that the non-interacting model (c^ = 0) is only compatible with WMAP 
data at 3a conhdence level, but when the supernovae data are included this is shifted to 2a 
conhdence level. 

We have shown that the Bayesian Information Criteria, that strongly disfavors increasing 
the parameter space to describe data sets of iV ~ 10^ points, provides positive evidence in 
favor of the existence of interaction. Other IQ models have been proposed Q] and even if 
their results are not directly comparable to ours, two interacting models, constructed with 
different motivations, suggest a value of the dark energy density smaller than in the concor¬ 
dance model, signals a need to investigate this type of models. As we have discussed, the 
quality of the £t and the cosmological parameters that can be derived by htting cosmological 
models to observations depend on the parameter space explored. 

The fact that IQ models appear to be favored by current observations suggests that dark 
energy and dark matter might not be so different entities after all. This is in line with 
recent ideas involving the Chaplygin gas. There, a single component plays the dual role 
of cold dark matter (at early times) and vacuum energy (at late times) and it interpolates 


between the two as expansion proceeds [^. The matter power spectrum could also be an 
important test of IQ models. A preliminary study shows that for the range of parameter 
space compatible with WMAP, the effect of the scalar held decaying into CDM has little 
effect 1^ . 

To summarize, the interacting cosmology model hts the WMAP data signihcantly better 
than the ACDM model does, and in fact alleviates the ISW effect at large angular scales, has 
no coincidence problem and provides a unihed picture of dark matter and dark energy. It 
predicts lower values of Hubble constant, dark energy density and higher baryon fraction. It 
is to be expected that the next generation of CMB experiments 0 and large scale surveys 
will enable us to constrain even further and discriminate between the different variants. 
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